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Abstract. This is a review of a decade-long effort to develop novel tools for explor- 
ing magnetism in cold astrophysical media and to establish a new field of molecular 
spectropolarimetry since Berdyugina et al. (2000). It is most directly applicable to the 
Sun, cool stars, substellar objects, planets and other minor bodies as well as interstellar 
and circumstellar matter. It is close to being a mature field with developed theoretical 
tools poised to uncover new insights into the magnetic structures in cooler environ- 
ments. Here I attempt a broad description of the literature and present some recent 
exciting results. In particular, following my programmatic review at SPW3, I discuss 
advances in molecular magnetic diagnostics which are based on the modeling of about 
a dozen diatomic molecules with various electronic transitions and states, including the 
most challenging - FeH. The applications stretch from sunspots to starspots, small- 
scale and turbulent solar magnetic fields, red and white dwarfs, and spin-offs such as 
polarimetry of protoplanetary disks and exoplanets. 



1. Advances in molecular diagnostics 



Interest in both the theory and experimental applications of the molecular Zeeman effect 
has had several ups and downs during the last hundred years. The peaks were around 
the 1930s, when the quantum mechanics was being established (e.g., |Kronig|1928[|Hill 



1929[ |Crawford| [T934), in the 1960s, when more advanced experimental facilities be- 



came available, and after the year 2000, when we discovered this niche for solar and 
stellar astrophysics (Berdyugi na et al.|[2000[ ). Since then we have learned to reliably 
model polarization of about a dozen diatomic molecules with various electronic transi- 
tions and states which are listed in Table [T] together with corresponding references. The 
upper and lower parts of the Table contain transitions which are either in the Zeeman 
or Paschen-Back regimes, respectively, at kG fields or lower. The references are also 
divided into two parts: the first lists theoretical work relevant for a given transition and 
the second lists the papers with modeling of particular lines under various conditions. 

A multitude of possible applications of molecular spectropolarimetry was first re- 
viewed at the SPW3 ( |Berdyugina et al.|2003b| ) and included solar and stellar magnetic 
fields at various scales. In this paper I have tried to collect a rather complete list of 
most relevant references for the past decade and review the recent progress in these 
applications and also present spin-offs from these developments such as polarimetry of 
protoplanetary disks and exoplanets. 
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Table 1 . Diatomic molecules having transitions in the optical and near IR that are 
present in solar and stellar spectra and are useful for measuring magnetic fields. 
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References: 1: Berdyugi na et al.| ( [2000| ), 2: |Berdyugina et al. (2001), 3: Berdyugina & Solanki 



(2001 



(2005 



2002 



Berdyugina & Solanki (2002), 5: Berdyugina et al. ('20033, 6: LBerdyugina et al. 
, 7: [Berdyugin a & Livin gston (2002), 8: Berdyugi na et al.| ( [2006a| ), 9: [Berdyugina et al. 
, 10: Berdyug inT&Flu ri (2004), 11: |Shapiro et al.]p007b| ), 12: |Shapiro et al.H2007a| ), 



Uitenbroek 



13: |Aframet al.| ( [2007l ), 14: |Afram et al.| ( [2008] ), 15: [Berdyugin a"etaLl ( [2007] ) 16: _ 
|et al.| ( |2004| ), 17: |Asensio Ramos et al.| ( |2004| ), 18: |Asensio Ramos et al.| ( |2005| ), 19: |Asensio 



Ramos & T rujillo Bueno (2006), +: see more references in Sect.[4| 0: this paper. 



2. 3D structure of sunspots 

Knowing the internal structure of sunspots allows us to gain insight on the energy trans- 
port in strong magnetic fields and, thus, on the processes inside the convection zone, 
where solar magnetic fields are generated and amplified before emerging at the sur- 
face on various scales. Spectropolarimetric studies of sunspots is a large field in solar 
physics, and there is already an immense amount of information collected during sev- 
eral decades (e.g. Solanki 2003 ). Nevertheless, there are still only a few studies on the 



3D structure of sunspots which involve depth-dependent inversions of Stokes profiles of 
multiple lines. The situation becomes even more complicated in spectra of the coolest 
parts of the umbra where molecular lines become too prominent to be ignored. This 
brings however the advantage to investigate umbra with a high sensitivity. For instance, 
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it was shown by |Mathew et al.| ( [2003| ) that simultaneous inversions of the infrared OH 
and Fe i lines at 1.56yum extends the sunspot model to higher atmosphere layers and 
greatly constrains temperature, magnetic field, and flows. Simultaneous inversions of 
the Fe i 5250A, TiO, and MgH lines holds promise to further improve average models 
of the sunspot umbra ( [Afram et al.|2006[|Arnaud et al.|2007t|Wenzel et al.|2010| ). 





Figure 2. Depth dependence of the temperature and magnetic field for different 
substructures of the sunspot umbra shown in Fig.[T] PD: penumbral dots, LB: light 
bridge, PUD: periferal umbral dots, CUD: central umbral dots, DU: dark umbra, 
SDU: darkest umbra. 
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It was also found that accounting for contributions of CaH and TiO lines in in- 
versions of the Fe i 6302A region is important for obtaining realistic models of the 
sunspot umbra (in prep.). Figure [T] shows our recent results obtained for the umbra of 
the sunspot NOAA 923 observed by Hinode from simultaneous inversions of atomic 
and molecular lines in the 6302A region. Thanks to the high spatial resolution of the 
the Hinode data, this provides a very detailed 3D structure of the umbra and allows 
us to extract average models of various substructures in the umbra (Fig. [2]). This can 
be compared with the MHD models of sunspots by (Remp eTet al.||2009| ) and provide 
further constraints for the simulations. Furthermore, realistic model atmospheres of 
sunspots are important for modeling solar irradiance variations which are modulated 
by sunspots on short-time scales. 



3. Imaging solar magnetic fields in molecular bands 



Imaging the solar surface in the G-band, where the main absorption is due to CH transi- 
tions, became a standard observing tool for studying the finest structure of small-scale 
magnetic fields and their dynamics (e.g., Berge r&"Title|200 j) . It was identified that the 
dissociation process of the CH molecule in hotter magnetic environments is the reason 
for their high contrast in such image s ([Steiner et al.|2001[ [Sanchez Almeida et al.|2001 



Schiissle r et al.||2003[ |Shelyag et al.||2004 ). Subsequently, we have shown that similar 



H: 



Ow 



Figure 3. Solar photosphere imaged in the U V OH (top) and CN (bo ttom) bands 
with the balloon-borne telescope Sunrise. From (Hirzberger et al.|2010j ). 



effects can be observed with other molecular bands, for instance in the violet systems 
of CN and OH (Berdyugi na et al.||2003a| ). Their dissociation equilibrium is similar to 
that of CH and one should expect a similar effect of a higher contrast in their bands. 
Indeed, ground-based imaging observations in the CN band at 3880A have confirmed 
our prediction that the contrast is even higher in the CN band as compared to the CH 
band ( [Zakharov et al.|2005[|Uitenbroek & Tritschler||2006[[2007] ). More recently, first 
imaging in the UV OH band were obtained with the balloon-borne Sunrise telescope 
(Hirzberg er et al.|[2010[ ). Theoretical interpretation and comparison with other bands 
will follow. 
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Figure 4. The NO A A 10667 sunspot imaged in the TiO band by V. Zakharov 
during our coordinated observations with THEMIS and SST in September 2004. The 
umbra region is shown with an enhanced contrast to reveal an amazingly complex 
network of umbral dots. 



Molecular bands which are stronger at lower temperature are useful for imaging 
temperature inhomogeneities in sunspots. For instance, it was predicted (Berdyugina 



et al.||2003a| ) that umbral dots should be seen at higher contrast if observed in the TiO 
7054A or UV OH band heads. Unfortunately there were no sunspots on the Sun during 
the Sunrise flight to verify this prediction, but ground-based observations at the SST, La 



Palma, in the TiO band were successful ( |Berger & B erdyugi na|2003| ). Figure|4]presents 



a sunspot imaged in the TiO band during our coordinated observing campaign with the 
SST and THEMIS (in prep.). It reveals a spectacular, never seen before, complexity in 
the umbra which is almost completely filled with umbral dots. This "web" of umbral 
dots appears to be a continuation of the penumbral structure but is deeply embedded 
in the umbra. This is in contrast to recent MHD simulations of sunspots by |Rempel 



et al.| ( [2009| ) where umbral dots are mere remnants of inhibited convection. Therefore, 



imaging and modeling of sunspots in the TiO band are important for understanding 
magneto-convection processes in the umbra. Furthermore, it is worth to mention that 
observing with the TiO filter apparently provides very sharp images of granulation and 
small-scale magnetic elements (but with a lower contrast than in G-band) and was even 
employed for recording Venus transit at SST ( http : //vt-2004 . solarphy sic s .kva. se| ) . This 
filter was also successfully employed for first-light, highest resolution images at the 



new 1.7 m telescope at the Big Bear Observatory ( |Goode et al.|2010| ). 



4. Turbulent magnetic fields on the Sun 

A study of turbulent magnetic fields on the Sun through scattering polarization and the 
Hanle effect is one of the major topics of the SPWs, and there are a number of contri- 
butions in this volume which address this topic. However, as they are devoted to atomic 
processes, I will try to cover here the key results obtained recently through molecular 
scattering polarization. The persistent questions concerning turbulent magnetic fields 
are (i) how strong are they on average, (ii) how are they distributed on the Sun, and 
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(iii) whether they vary with time, e.g., with the solar cycle. It appears that the differ- 
ential Hanle effect applied to molecular lines has the greatest potential to answer these 
questions. 

The first clues to puzzling molecular polarization in the second solar spectrum 
were provided by Ber dyugina et al.| ( |2002| ). Even though in that work we did not discuss 
the effect of upper level lifetimes on the sensitivity to the Hanle effect, as was pointed 
out by |Trujillo Bueno| ( [2003] ) and |Landi Degrinnocentil ( [2003|[2007] ), it a ppeared that 
the majority of lines within the same molecular band have similar lifetimes ([Berdyugina 



& Fluri 2004), so this argument does not really clarify the question of why molecular 



polarization appeared nonvariable. 
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Figure 5. Scattering polarization in two C2 triplets. Two fits (solid and dash- 
dotted) to a synoptic observation (dotted) (averaged along t he slit) are obtained using 
the differential Hanle model by Berdyugina & Fluri ( 2004). The profiles shown with 
the solid line were calculated using a magnetic field of 7 G and those shown with 
the dash-dotted line represent 90 G. There is a clear indication that the field cannot 
exceed a few tens of Gauss but its value can only be determ ined if one assume s that 
the two triplets are differently affected by collisions. From|Kleint et al.|([2010b). 



Answering the first question above is a nontrivial task. The strength of the entan- 
gled magnetic field strongly depends on its assumed distribution within the resolution 
element on the Sun as well as on the model atmosphere used (e.g., |Trujillo Bueno et al 



2004| ). Most results so far were obtained assuming an isotropic distribution, which 



provides a simple analytical tool to infer the field strength (Stenflo 1982). When this 
is applied to lines with similar thermodynamical properties but with different mag- 
netic sensitivities, the differential Hanle effect provides the most model independent 
field strength. At the SPW3 |Trujillo Bueno| ( |2003| ) argued that the observed differ- 
ence in scattering amplitudes within partly resolved triplets of C2 transitions (such as 
the P-triplet at 5 141 A in Fig. [5} in the second solar spectrum are due to the Hanle ef- 
fect. However, it could also be explained by a simple blending effect of two unresolved 
triplet components as compared to a partly resolved one. In search for the truth, we have 
carefully investigated the atlas of Gandorfer ( 20 00]) and found a fu lly resolved and prac- 
tically unblended R-triplet at 5 140 A ( |Berdyugina & Fluri||2004[ see also Fig. [5]). This 
triplet was clearly affected by the Hanle effect and not the way it was argued before. 
However, even in this clear case the interpretation was hindered by unknown parame- 
ters, such as Einstein coefficients, which are not trivial in the case of C2 transitions, and 
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Figure 6. L<?/h 2/7 amplitude ratios for 78 synoptic measurements of the C2 
triplets shown in Fig. [5] sorted by heliographic position angle. The observed varia- 
tions are within one standard deviation for all angles. This indicates a surprisingly 
homogeneous distribution of turbulent magnetic fields during the solar minimum of 
2008-2009. Right: Amplitude ratio variations with time. A systematic decrease of 
R2/P and increase of R1/R2 and R3/R2 seem to occur as the solar activity declines. 
Second-order polynomial fits to the line ratios are shown with dotted lines. From 
|Kleintetal.| ( |2010b| ). 



collision rates, which are not known with the necessary accuracy (e.g., |Asensio Ramos 
& Trujillo Bllenol|2005t |Bommier et"aT1[2006t |Kleint et aL]|2010b| ). These issues are 
not still clear and need further investigation. However, to summarize, it appears that 
many studies of the Hanle effect in molecular lines deduce a relatively low turbulent 
magnetic field strengths o f 5-30 G ([Faurobert et al.||2001 ; Faurobert & Arnaud|| 2003 ; 
Berdyugina & Fluri|[2004HDerouich et aL||2006[|Trujillo Bueno et al.|2006[ |Bommier 



et aL|2006[|Shapiro et al.|2007a[|Kleint et aL|2010b| ), which may be an indication for a 
depth-dependent effect. It is clear however, as illustrated in Fig. [5} that the field strength 
in the lower atmosphere cannot exceed a few tens of Gauss, in contrast to conclusions 
by, e.g., |Trujillo Bueno et al.1 ( [2004| ). 

In such a situation the second and the third questions cannot be answered without a 
systematic and homogeneous study. That is why we have initiated a synoptic program at 
IRSOL for measuring during several years (hopefully through at least one solar cycle) 
the scattering polarization in the C2 triplet and a few atomic lines. The results for 
the first two years covering the solar minimum in 2008-2009 were reported by Kleint 
et al.| ( [2010b[ |201 1] ). Figure [6] illustrates that turbulent fields were evenly distributed 
over the Sun during this solar minimum. The field strength was on average at least 
4.7+0.2 G (without collisions taken into account) or 7.41+0.09 G (with collisions and 
a NLTE analysis). It did not vary during the last two years, but there could be a very 
small variation of the turbulent field strength (3<x-limit) since the solar maximum in 
2000 (based on data of Gandorfer 2000). This result is however very sensitive to the 
available statistics, so we need to continue the synoptic program to clarify whether 
turbulent magnetic fields are generated by a global dynamo responsible for the solar 
cycle or a local dynamo which may not vary at all on the cycle time scale. It is worth to 
mention that such a program would be a rather minor expense for large solar telescopes 
such as ATST (see |Kleint et aL|2010a| ). 
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5. 3D Structure of Starspots 



Sunspots are the primary evidence of solar activity, and they harbour the strongest 
field and the coldest plasma on the solar surface. Such a phenomenon is common 
among stars possessing convective envelopes, where magnetic fields are believed to 
be generated ( [Berdyug ina 2005). However, magnetic fields were never measured in- 
side starspots, and their internal structure remained unknown (e.g., |Berdyugina||2009[ 
Strassmeier|2009[|Donati & Landstreet|2009| ). Studying this phenomenon is especially 



interesting for cool dwarfs, which are known to show strong magnetic activity and 
whose dynamo is apparently different from that of more massive cool stars, such as 
the Sun. The process of generation of strong magnetic fields in objects with very deep 
convection zones or even fully convective is not yet understood. We have developed a 
unique approach to probe their atmospheres and magnetospheres with polarimetry and 
reveal a 3D structure of starspots. 

The prospects of molecular spectropolarimetry for exploring magnetic interiors 
of starspots were foreseen by Berdyugina (2002), but adequate observational facilities 
were not available at that time. Once the high-resolution spectropolarimeter ESPaDOnS 
( [Donati et al.|2006| ) became available at the CFHT, our team obtained first polarimetric 
measurements of starspots on M and K dwarfs. We have detected circular polariza- 
tion in many lines, including the first detection of polarization in TiO, CaH, and FeH 
transitions of 0.5-1% with the noise level of 0.1% (Berdyugina et al. 2006b) and in 
chromospheric emission lines ( Berdyugina et al.||2008b ). Here we summarize first re- 
sults inferred from these data (in prep.) for an active red dwarf AU Mic, MlVe, which 
is the brightest X-ray source among such stars. 

We selected several lines of Fe i and Ti i and many of TiO, CaH and FeH with 
prominent Stokes V profiles (Figs.|7]&[8]). These lines cover a wide range of tempera- 
ture and magnetic sensitivity and also effectively form at different heights and locations 
in the stellar atmosphere. This provides a unique opportunity for probing directly the 
interior of starspots, basically independent on the distance to the star or its size, i.e. 
clearly beyond current direct spatial resolution on stellar surfaces. Since the observa- 
tions of AU Mic were obtained during two consecutive nights, which roughly covers 
half the stellar rotation period, we can only analyse spots visible during that time. We 



applied our new Zeeman-Doppler imaging (ZDI) inversion technique (Sennh auser et al 



2009b) and Stokes parameter synthesis ( [Berdyugina et al.|2003a| ) to individual species 



forming at different heights in the range of 60-2 10 km above the photosphere. Due 
to lack of spatial information in the line profiles (because of the slow stellar rotation, 
vsin/=10 km/s), shapes of spots remained uncertain and for simplicity were initially 
assumed to be circular and homogeneous, which was subsequently refined by inver- 
sions until the best fit to the data was achieved. The five obtained surface images are 
shown in Fig. [9j 

Spots on AU Mic were found in a near-equatorial zone (their exact latitude is 
uncertain because of the unknown inclination angle of the stellar rotation axis). This 
is similar to the Sun but in contrast to rapidly-rotating active stars where Coriolis force 
brings spots closer to the rotational axis, i.e. poles. Even though the spot shapes are 
unknown, their area and longitudes are well constrained by the data. We identify four 
spots, one of positive and three of negative polarity (Fig. [9]). They are 500-700 K cooler 
than the photosphere and harbor a maximum magnetic field of 5.3 kG. The smallest spot 
is comparable in absolute dimensions (70 Mm) with very large sunspots (60 Mm), and 
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Figure 7. Stokes profiles of Fe i (two left panels)) and Ti i (three right panels) 
lines in the spectra of AU Mic. Data (crosses) for four rotational phases recorded 
during two consecutive nights are shifted in vertical for clarity. Model fits (red line) 
correspond to the images shown in Fig. [9] 
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Figure 8. The same as Fig. [7] for molecular lines. 
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Figure 9. Starspots on AU Mic at different heights above the photosphere obtained 
from different lines. Darker regions are cooler (scaling is individual for each image). 
Green and orange sticks depict magnetic vectors of away and towards the observer, 
respectively. They are scaled to the unit length corresponding to 5 kG. 



it practically disappears at the bottom of the atmosphere. The vertical temperature and 



magnetic field distributions inside one starspot are shown in Fig. 10 



The spot temperature decreases with height from 3100K to 2400 K (+100 K), 
with the gradient dT/dh~5 Kkm" 1 (Fig. 10), which is significantly larger than that in 
sunspots. The magnetic field strength (modulus) in spots reaches more than 5 kG, which 
is much stronger than the equipartition level. Its radial component (which was found 
to be dominant in spots) decreases with height from 4.0 kG to 2.3 kG (+0.5 kG). The 
gradient <iB r /d/z~10Gkm _1 is an order of the magnitude larger than that in sunspots. 
To exclude the influence of the density scale height, we consider the gradient dB r /dT, 
which is about 2GK" 1 in AU Mic spots and about 1 GK" 1 in sunspots, and they are 
both significantly larger than the adiabatic gradient. It appears that this difference oc- 
curs due to decoupling of the magnetic field from the ambient gas, whose ionization 
degree is much lower on AU Mic than on the Sun. This goes well along with the fact 
that the spot area on AU Mic increases with height, i.e. the field fans out and is loosely 
confined by the gas. 

We find also evidence for the second magnetic component in the AU Mic atmo- 
sphere: a mix-polarity field of 3.5+0.5 kG covering about half of the photosphere out- 
side cool spots. This field is apparently confined to the lower photosphere or perhaps 
weakens with height. It appears analogous to a solar network magnetic field. Its topol- 
ogy is unknown and for simplicity assumed to be radial. This result clarifies high mag- 
netic filling factors deduced from measurements of Zeeman-broadened profiles, includ- 
ing FeH, in active M dwarfs ( |Johns-Krull & Valenti|1996 , Valenti et al.||2001[ Reiners 



& Bas ri 2006 ). With our advanced quantum-mechanical modeling of FeH transitions, 



the accuracy and reliability of such measurements has significantly improved (Afram 
|etal.|2009l|20TTl). 



Less active targets show however polarization signals significantly smaller than 
1%, while current observational facilities limit the accuracy per wavelength to 0.1%. 
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Figure 10. 3D structure of a starspot on AU Mic (vertical cut through the atmo- 
sphere). The temperature and radial magnetic field decrease with height, while the 
spot area increases. An entangled field of 3.5 kG (assumed to be radial) fills space 
between spots in the lower photosphere. 



Therefore, for such targets it is not possible to detect signals in individual lines, and 
multi-line techniques such as LSD ( |Donati et al.|1997| ) are commonly used to increase 
the signal-to-noise ratio (SNR). However, limitations of the LSD are now well recog- 
nized (e.g., |Berdyugina||2009| ), and new techniques emerge to advance our ability to 
measure and interpret Stokes profiles, such as PC A (Martinez Gonzale z et al.| |2008 ), 
NDD ( [Sennhauser et al.|2009a| ), and ZCD ( [Sennhauser & Berdyugina|2010b| ). The lat- 
ter two are particularly powerful as they take into account nonlinear effects in blending 
spectral lines and their Zeeman components for both atomic and molecular transitions. 
At the moment the ZCD is the most sensitive technique for detecting very weak mag- 



netic fields, such as about 0.5 G on the K giant Arcturus ( [Sennhauser & Berdyugina 
2010a| ). When combined with the ZDI, it will open new opportunities to explore mag- 
netic spots in many objects across the HR-diagram. 



6. White Dwarfs 



Applications of molecular spectropolarimetry extend even to cool white dwarfs (WDs), 
where this technique is the only hope to clarify the frequency of magnetic objects and 
learn about their evolution. For instance, a small group of cool, helium-rich WDs 
(DQ type, T<7000 K) show only carbon-based molecular bands, mainly C2 bands (e.g., 
Dufour et aL||2005| ), which are therefore the only indicators of their magnetic fields. 
It is believed that the majority of these stars are not magnetic, since no polarization 
was reported in C2 until now. However, a famous representative of this group, G99-37, 
for more than four decades having been a unique WD with both CH and C2 molecular 
bands, is clearly magnetic, as indicated by prominent broad-band polarization in the 
CH bands discovered by Ange l & Landstreet| ( |1974[ ). This is due to magnetic dichroism 
appearing in molecular bands in the presence of a strong magnetic field, which perturbs 
the internal structure of the molecule and results in net polarization due to the Paschen- 
Back effect (PBE). Recently we analyzed new spectropolarimetric observations of G99- 
37 ([Berdyugina et al.|2007| ) using our approach to the molecular PBE ( |Berdyugina et al 
2005| ). In addition to previously known molecular bands of the C2 Swan and CH A - X 
systems we found also a firm evidence for the violet CH B - X bands at 390 nm and 
C2 Deslandres-d'Azambuja bands at 360 nm. Combining the polarimetric observations 
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Figure 1 1 . Observed intensity distribution, circular and linear polarization of the 
white dwarf GJ 841B - the second WD with CH bands (thin line). The CH and C 2 
bands are identified. Synthetic spectra for a longitudinal field of 1.3 MG directed to 
the observer shown by thick red lines. 



with our model calculations, we deduced a dipole magnetic field of 7.5+0.5 MG with 
the positive pole pointing towards the Earth. 

Inspired by this result, we started a spectropolarimetric survey of cool DQ-type 
WDs. This has recently resulted in discovery of the second WD with polarized CH 
bands, in both A - X and B - X systems! The CH bands in the intensity spectrum 
of GJ 84 IB are completely masked by strong C2 absorption, and only polarimetric 
measurements of high accuracy could reveal the CH bands and the magnetic nature of 
the WD ( Vornane n et al.|2010| ). By modeling four Stokes profiles in the CH bands we 
determined the magnetic field strength of 1.3+0.5 MG, which is largely longitudinal 



(Fig. 11 ). It is interesting that the circular polarization is of the same sign as in G99- 



37, and the continuum polarization of +0.2% is also present. In addition, GJ 841B is 
the first white dwarf to show signatures of weak polarization in the C2 Swan bands. 
The discovery of the second WD with polarized CH and C2 bands holds promise that 
there are more such objects out there. The developed technique of molecular magnetic 
dichroism has proved itself to be an excellent tool for studying magnetic fields on cool 
and faint magnetic objects. 



7. Exoplanets and Protoplanetary Disks 

Two interesting spin-offs from our developments for scattering polarization in the solar 
atmosphere are polarimetric studies of exoplanets and of protoplanetary disks. Here 
we summarize some of the major results obtained for exoplanets in recent years. The 
development of a new diagnostic technique based on radiative pumping of absorbers to 
explain ubiquitous absorption line polarization in embedded stars and to explore inner 
parts of protoplanetary disks has been pioneered by |Kuhn et al.| ( |2007] ) and is reviewed 
in these proceedings by |Kuhn et al.| ( |20lT] ). 
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Figure 12. Left: A star emits unpolarized light in all directions. The part of light 
that is scattered by a planet towards the observer is linearly polarized perpendicular 
to the scattering plane. Thus, the observer registers a partially polarized light from 
an unresolved star-planet system and can detect directly the scattered light using 
a polarimeter. Right: As the planet revolves around the star, the scattering angle 
changes and both the direction and degree (red sticks) of polarization vary. For 
a circular orbit, as shown here for HD 189733b, two peaks per orbital period can 
be observed near elongations. Therefore, polarimetry provides information on the 
planetary atmosphere outside transits. 



The physics of scattering in the planetary atmosphere is similar to that of the solar 
case. The light scattered off the planet is linearly polarized perpendicular to the scatter- 
ing plane (Fig.[l2j left). In general, when the planet revolves around the parent star, the 
scattering angle changes and the Stokes parameters vary. If the orbit is close to circular, 
two peaks per orbital period can be observed (Fig. |T2j right). Thus, the observed po- 
larization variability exhibits the orbital period of the planet and reveals the inclination, 
eccentricity, and orientation of the orbit as well as the nature of scattering particles in 
the planetary atmosphere. It is feasible therefore to detect the scattered light with po- 
larimetry, even if an exoplanet is not spatially resolved, and hot Jupiters on close orbits 
are the best candidates for polarimetric detections. 

Our first target was a very hot Jupiter HD189733b which is only 0.03 AU away 
from the star. In 2006-2007 we observed the system in the 2?-band with the double 
image CCD polarimeter DiPol installed at the remotely controlled 60 cm KVA telescope 
on La Palma ( |Berdyugina et al.|2008a| ) and in 2008 in the UBV-bmds with the TurPol 
polari meter at the 2.5 m Nordic Optical Telescope (NOT), La Palma (Berdyugin a et al 



2010). Our observations clearly reveal polarization peaks of 10 4 near elongations as 



expected for a planet on a circular orbit (Fig. [13]). 

To interpret the data, we constructed a semi-empirical model atmosphere where 
polarization is caused by Rayleigh scattering ( |Fluri & B erdy ugina [ |20 1 0[ ) . From the 
phase-locked polarization variations we determined orbital parameters, including the 
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Figure 13. Polarimetric data (normalized Stokes q and u) for HD 189733b in U 
(pink) and B (orange and purple for two different sets of data) photometric bands 
(upper panels). The binned data are shown in lower panels (blue). The mean error 
of the binned data is 1.7-10" 5 . Curves are the best-fit solutions for a model atmo- 
sphere with Rayleigh scattering on dust particles. The standard deviation of the fit is 
1.110" 5 . 



positional angle on the sky plane (see Fig. [12} right). Our model atmosphere includes 
the following opacity sources: (i) scattering on H, H2, He, CO, H2O, CH4, electrons, 
MgSi03 dust condensates, and (ii) continuum absorption due to free-free and bound- 
free transitions of H, H", H+ , H~, He, He", Si, Mg, and Fe. To obtain the best fit 
to our data, we varied temperature, gas pressure, and dust particle size and density 
TVd in the planetary atmosphere within the height range of 82,000-84,500 km with a 
step of 100 km. For each layer we calculated opacities and polarization and searched 
for the best- fit solution with a x 2 minimization procedure. The dust distribution was 
considered as a layer with a Gaussian profile, which is characterized by the peak values 
of and In particular, we found (in prep.) that there is a set of similar quality 
best-fit solutions with a relation between the peaks of and (Fig. [T4j left). 

The data taken at different wavelengths including upper limits by others clearly 
indicates the dominance of Rayleigh scattering in the optical and can be described by a 
single model atmosphere (Fig.[l4| right). A transit curve calculated taking into account 
opacities in the atmosphere perfectly agrees with the transit data by |Pont et al.| ( |2007"] ). 
The geometrical albedo of our model atmosphere varies with wavelength and reaches 
0.48 at 400 nm, 0.2 at 550 nm, 0.1 at 600 nm, and 0.05 at 700 nm, in agreement with 
upper limits for hot Jupiters. We conclude that the polarization in HD189733b occurs 
due to Rayleigh scattering in the atmosphere, most probably in on dust condensates of 
average size of about 20 nm. 

Our results establish polarimetry as an important tool for studying directly exo- 
planetary atmospheres in the visible. In the near future it will be employed for non- 
transiting systems. 



8. Conclusions 

This review demonstrates that polarimetry of cool atmospheres with the help of molec- 
ular physics is a maturing field with many theoretical tools developed and ready to be 
employed to uncover new insights into magnetic structures and inhomegeneities of cool 
atmospheres. In addition, as we learn about the physics of the processes responsible for 
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Figure 14. Left: x 2 -contours for peak values of the dust particle size and num- 
ber density The confidence intervals of 68%, 90%, and 99% are shown with 
solid lines. The dashed red line indicates a relation for acceptable solutions. Right: 
Stokes q pe ak amplitudes in diffe rent passbands. Open circle s are the UBV ampli- 
tudes from (Berdyugin a et al.|2010| ), cross is the upper limit by Wiktorowicz ( 2009 ), 
and triangle is the upper limit for hot Jupiters by |Lucas et al.| ( |2009| ). Horizontal bars 
indicate the passbands' FWHI. Solid red line is the polarization curve for the Model 
D2 with max a d =50nm and max N$= 5-10 6 cm -3 . Dashed and dashed-dotted curves 
are the same model but with max of 40 nm and 60 nm, respectively. Our data 
are consistent with other measurements and indicate that polarization occurs due to 
Rayleigh scattering with a peak in the blue. 



polarization in molecular transitions, a great potential for interdisciplinary research is 
revealed. For instance, polarimetry of exoplanets and protoplanetary disks are new 
emerging fields where this expertise is of great advantage. 
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